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Abstract 

This work focuses on an unexplored aspect of non-symmetric geometry where only the off- 
diagonal metric components along the extra dimension, in a 5-dimensional spacetime, are non- 
symmetric. We show that the energy densities of the stationary non-symmetric models are similar 
to that of brane models thereby mimicking the thick-brane scenario. We find that the massive 
test particles are confined near the location of the brane for both growing and decaying warp 
factors. This feature is unique to the non-symmetric nature of our model. We have also studied 
the dynamical models where standard 4D FLRW brane is embedded. Our analysis shows that the 
non-symmetric terms deconfine energy density at the early universe while automatically confine at 
late times. 
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I. INTRODUCTION 



General relativity (GR) has survived numerous experimental tests [l| and achieved the 
status of the accepted theory of gravitation. Though observations, relating to galactic 
rotation curves 0, 3] and supernova data 4| can only be explained within the framework of 
GR with exotic matter fields. These lead to inception of many new ideas in cosmology as 
attempts to justify the observational facts. In the literature, broadly, there have been two 
approaches: (a) Assume GR is correct and include exotic matter fields in the right hand side 
of Einstein equation (like dark matter or dark energy) and (b) Include the possibility that 
GR is valid only upto a certain length scale. At largest and smallest length scales, theory 
of gravity is modified. Such models can potentially avoid introduction of exotic matter 
fields and can plausibly explain observations. Until now, both these approaches have been 
partially successful to explain observations. However, any modification to GR has to be such 
that it does not contradict with GR in the regime where it is consistent with observations. 

Most of the modifications to GR has focused on introducing higher order curvature terms 
which dominate at larger length scale or smaller length scale. However very few attempts 
have been made in the literature to look at the simplest of plausible extensions to the 
assumption of GR- non-symmetricity of the affine connection. It is important to understand 
that GR assumes that spacetime manifold is metrically connected. GR further assumes 
that the metric tensor is symmetric which eventually leads to symmetric affine connections. 
Physically this means that the corresponding spacetime is torsion free. However, geometric 
nature of gravity does not require that. There have been attempts to construct, more 



fundamental, affine gravity models in the literature 



In reference 



a 



6], Schrodinger described the field equations of non-symmetric purely affine 



gravity theory. Moffat has extensively studied the, so called non-symmetric gravitational 



theory (NGT) for over four decades 8|- 14[. One of the key results of NGT models, unlike 
GR or hig her derivative gravity models, the black- hole solutions are singularity free [s, si 
12L Il5j . Recently, attempts were made to show that NGT can be used as an alternative 



11 



to dark matter and dark energy 16Ml8j|. Attempts have been made to constraint Moffat's 



NGT through observations using the results of recent free-fall Galilean test of the WEP 



19l | . However this class of models has some difficulties related to ghost poles 20] and 



non-uniqueness [211 ]. 
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Models with extra dimensions (D > 4) have also been studied in detail [22] and gained 



substantial attention, in the last few decades, mostly motivated from string theory 



231. In 



the last decade or so, warped braneworld models are extensive 



y analysed, given its potential 



to solve some longstanding puzzles in theoretical physics [24j-|26|. Thick brane solutions (3- 
brane with, finite width, embedded in a 5D bulk) of Einstein equations in warped spacetimes, 
which are natural generalisation to original infinitely thin brane scenario of RS, appear in 
various multidimensional field theories coupled to gravity {27 1. Note that brane thickness is 
an essential ingredient for extending the idea of brane world to multidimensional spacetimes. 
One of the key assumptions in the brane models is that the 5 th dimension has a discrete 
symmetry which distinguishes them from the other 4 dimensions indirectly, explaining the 
non-observability of the extra dimension. While this is a crucial feature, it is not clear how 
to include this discrete symmetry from fundamental principles. 

In this work, we show that the continuous symmetry i.e. non-symmetric metric for 
the higher dimension mimics the feature of thick brane scenario. We analytically construct 
stationary and dynamic solutions of field equations in simple five dimensional non-symmetric 
bulk spacetimes, where the contracted torsion tensor vanishes, and analyse the corresponding 
energy densities and timelike geodesic motion. We assume that the (off-diagonal) metric 
components along the extra dimension is anti-symmetric while the metric components of 
the 4D part is symmetric. 

We find that the geometry itself confines the timelike geodesies near the location of 
the brane for both growing and decaying warp factors. This feature is different from the 
symmetric bulk models where confinement happens only in presence of growing warp factor 
. The dynamic solutions has the standard 4D FLRW universe embedded in the 5D bulk. 



This paper is organised as follows. In section II, we briefly review and introduce the field 
equations in non-symmetric gravity. In sections III and IV, we analyse the stationary and 
dynamic solutions respectively. Section V has a brief discussion over possibilities of finding 
dynamic warped braneworld solutions. Then we conclude with a summary of the results 
found in Section VI. Signature convention used is (—1,1,1,1,1). Capital Latin indices, 
Greek indices and small Latin indices are used to denote the coordinates of 5D spacetime, 
4D spacetime and 3D space respectively. Ordinary and covariant differentiations are denoted 
by a coma and a semi-colon, respectively. An overdot Q means differentiation w.r.t. the 
affine parameter A. 
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II. NON-SYMMETRIC GEOMETRY 



The gravitational field equations in non-symmetric spacetime, in (D + 1) dimensions, are 



given by 



6] 



-99 



[AB] ),b = 0, (2.1) 



9ab,o ~ fD AC 9os ~ r D OB 9 AD = (2.2) 

G AB = R AB ~ 2 9 *B R = 87rT AB ( 2 - 3 ) 

where g AB and T AB are the metric and matter tensors respectively. '[*]' means anti- 
symmetrisation ('(*)' will mean symmetrisation). r^ c 's are the modified affine connections. 
R and R AB are the Ricci scalar and the first Ricci tensor, respectively, which are defined as: 

R = g AB R AB where R AB = R C ACB (2.4) 

and the Riemann tensor is given by 

R D _fD _ pD , f K f D _ f K f D (n ^) 

1X1 ACB ~ 1 AB.C AC,B AB KC ^AC^KB- 

Note that, since the spacetime is non-symmetric, the ordering of indices in the definition of 
the tensors is important 6|. 

The second Ricci tensor is defined as 

Qab = r° CB , A - r° CA , B (2-6) 

which is anti-symmetric and has the form of a curl. The covariant derivative of a vector S A 
is defined as 

where T are the original affine connections of the manifold. Note that T c and r c are 

BC & [AB] [AB] 

both (torsion) tensors, however the connections are not. The modified affine connections, 
T A are recursively related to the affine connections T A as 

BC a BC 



f A =T A --5iT c where T c = -(T D — T D ) (2.8) 

BC BC ]J 2 



which implies 



rc = ^ D -f^) = 0. (2.9) 



Note that Eq. (I2.9p is identical to Eq. ( 12. ip . Since Eq. ( 12. 8 p is recursive, it is difficult to 
solve analytically. Hence, we assume the following 

= i(C - C) = 0, (2.io) 

which essentially implies 

This further implies that the second Ricci tensor, Eq. (12. 6p . also identically vanishes. Most 
importantly, with this assumption we reduce the solution space to one of its subclasses. We 
show that Eq. (12.111) leads to rich class of interesting solutions, specifically with only one 
space dimension being non-symmetric, following the steps below: 

1. Assume a simple 5D metric ansatz whose 4D part is symmetric and the off-diagonal 
components corresponding to the extra dimension is non-symmetric. 

2. Derive the affine connections using Eq. ( 12. 2ft . 

3. Then determine the metric functions solving the constraint equations (I2.9P or (I2.10p . 



4. Analyse the resulting matter energy density G o assuming that these models are solu- 
tions (numerically) of field equations with appropriate matter energy tensor added to 
the right hand side of Eq. (12.3ft . 

5. Solutions of the corresponding geodesic equations then reveals how motion of test 
particles is modified by the non-symmetric nature of the space time. Investigation of 
whether Weak Energy Condition, G AB u u B > 0, is satisfied or not, tell us to what 
extent these models are compatible with general relativistic framework. 

In the next Section, we discuss stationary spacetimes. 

III. STATIONARY BULK SPACETIME 

Let us assume the following simple ansatz for the stationary 5D metric: 
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/-/Wo o 

O f(a) 

f(a) 





d(a) \ 
d{a) 
d(a) 
f(a) d{a) 

) 



(3.1) 



\-d{a) -d{a) -d{a) -d{a) r' 

where the non-zero metric components are functions of the extra spatial dimension a. Here 
r is a constant and in all the figures, presented later, it's value is taken to be unity. Note 
that the line element corresponding to metric (13. ip . 



ds 2 = f(a)[-dr] 2 + dx 2 ] + r 2 da 2 



(3.2) 



does not capture the true nature of the bulk which is stationary. 

Solving Eq. (12. 2p with our ansatz, the components of can be obtained as functions 
of f(a) and d(a) (see Appendix). The constraints (I2.10p imply 



d 



3^/, CT + rM^ + ^ 



0. 



(3.3) 



P " ' \f ' d 

The above equation has two unknown metric functions, hence, one has to choose one and 
determine the other. Note that the results found does not depend on the sign of d(a) since 
one can always exchange the diagonally opposite anti-symmetric metric components. 
We will study two different scenarios - A) f(a) = cosh 2 (cr) and B) f(a) = sech 2 (cr). 



A. f(a) = coshV) 

Putting f(a) = cosh 2 (cr) in Eq. ( 13.31) one gets: 

d(a) = ± r (3.4) 

^dr 2 coshV) -2sechV) 

where C\ is the integration constant. Note that d(a) can be real or purely imaginary, for 
all a, depending on whether C\ > 2/r 2 or C\ < respectively. For < C\ < 2/r 2 , d(a) is 
singular at \a\ < cosh~ 1 (2/C7ir 2 )( 1 / 6 \ Fig. CQdepicts energy density profiles in two such cases 
where d(a) is either real (r = 1, C\ = 3) or imaginary (r = 1, C\ — —1). Energy densities 
in both the situations behave almost similarly i.e. reaches maximum at o = 0. In this 
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sense, these models mimic a thick brane located at a = 0. This behavior is not drastically 
different from the symmetric case because with increasing |cr|, d(o~) becomes negligibly small 
compared to f(o~). Note that this case mimics the thick brane scenario, with a growing warp 




a- a cr 



FIG. 1: Plots of energy densities vs a for growing warp factor with d{a) real, imaginary and zero 
respectively. 



factor, where the brane is located at a = 0. 
B. f(a) = sechV) 

With f(cr) = sechV), Eq. implies: 

d(a) = ± r (3.5) 

^C 2 r 2 sech 4 ((x) -2cosh 2 ((j) 

where C 2 is an integration constant. Note that d(a) is purely imaginary, for all a, when 
G\ < 2/r 2 . For C\ > 2/r 2 , d(a) is real below \a\ = cosh~ 1 (2/Cir 2 )^ 1 / 6 ) where it becomes 
singular. Fig. [2] depicts energy density profiles in two such cases where d(a) has singularities 
(r = 1, C*2 = 5) and purely imaginary (r = 1, Ci = — 1). It is interesting to note that non- 
symmetric components does drastically change the energy density profile compared to the 
symmetric scenario, unlike the previous model where the effect of d(a) is negligible for large 
a. This case mimics the thick brane scenario, with a decaying warp factor, where the brane 
is located at o = 0. Note that when d(a) is purely imaginary, energy density has a local 
maxima at a = and the location of the brane is less stable compared to the other situations. 

In the next section, we look at the corresponding geodesic equations to understand the 
nature of the spacetime discussed in this section 
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FIG. 2: Plots of energy densities vs a for growing warp factor with d(a) real, imaginary and zero 
respectively. 



C. Geodesies 



The geodesic equations, satisfying the time-like constraint 

/(-7) 2 + ^^) + r 2 d 2 = -l, (3.6) 

i 

derived for general ansatz can be written in the following simplified form: 

V + - YJT^ + x i) & = °> ( 3 - 7 ) 

J i 
Xi + ^j-Xi& + (?) + ^Xj)cr = 0, (3.8) 

J i 

where F = (2d 2 j f + r 2 ). Eqs. fl3.7l) - fl3.9p imply that non-symmetric components do effect 
the geodesic motion although they do not appear in the line element (13. 2\i . Note that as 
r] and Xi are the cyclic coordinates, one expects f] oc 1// and ii oc 1/f. However such 
solutions does not exist except in special cases when the last terms in the left hand side of 
the geodesic equations vanish. 

This could occur for three cases- 

1. / oc d 2 : this condition along with Eq. (13 .3p would imply f,d = constant. 

2. f] = Xi = 0: implying there is no flow of time which would be unphysical. 

3. T] = — Ylii x i'- this corresponds to set of geodesies in symmetric spacetime. In general, 
putting F = constant scenario refers to the geodesies in a symmetric 5D bulk, that 



are studied in detail in 



28] 



Note that in our case behaviour of geodesies can not be understood from the so called 



geodesic potential 



29| derivable from the line element as such. As already pointed earlier, 



effects of non-symmetric nature of the spacetime does not show up in geodesic potential. 
Here we are interested in the role played by the non-symmetric components of the metric. 
The equations (I3.7ti3.9p are highly coupled and hence difficult to solve analytically. Hence 
we solve them numerically and plot the results. 

Note that the qualitative aspects of the geodesies do not change with the variation of 
initial conditions. 

1. f((x) = coshV) 



The geodesies in this growing warp factor scenario are shown in Fig. [3j In all the cases 
geodesies look similar and essentially their motion along the extra dimension are confined 
to a = 0. Effect of the non-symmetric nature do not modify the cx-component. However 
when b(a) is imaginary £ r components also become oscillatory. 



f(u)=coslr(tr). d(<r)=±|3cosh V)-2sech 2 ([j)] 



x(A) 



hin^a ls ]|-un. dlin— ilcush !ir)+2seclr[ir)| 



x(A) 



f(<7-)=CoshV), d( t r)=0 



x(A) 

! —a-m 



riOi) 



FIG. 3: Plot of geodesies Xj(A) and cr(A) vs 77(A) for growing warp factor with d(a) real, imaginary 
and zero respectively. Initial conditions are Xi(0) = <r(0) = 0, r/(0) = and ij(0) = (i(0) = 1, r)(0) = 
y/E. 



The oscillatory behavior of cr(A) with d(a) = 0, can be analytically understood as follows. 
With growing warp factor Eq. f)3.9p . to the leading order in a can be written as 

a + axC = (3.10) 

where Cq is a positive quantity. Eq. (I3.10P clearly resembles equation of simple harmonic 
motion. However this equation is non-trivial for d(a) 7^ 0. Numerical analysis shows that 
the new coefficient appearing in the geodesic equations also contributes in a similar manner 
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in this case. From Fig. HJ the oscillation of cr(A) for real d(o~) seems obvious. Whereas 
for imaginary d(a), the oscillatory behaviour of cr(A) stems through £j(A) as in Eq. ( 13. 8p . 
the non-symmetric contribution (or F 2 in Fig. HJ), comes with a negative sign. It is worth 
mentioning here that these reasonings also establish the reliability of the numerical results 
presented. 



f(o-)=cosh 2 (cx), d(o-)=+[3cosh 4 (o-)-2sech 2 (o-)J 2 




f(o-)=cosh 2 (cr), d(o-)=+i|cosh 4 (o-)+2sech 2 (cr)] 



F 2 



-3 -2 



2 3 



FIG. 4: Comparing the coefficients appearing in the geodesic equation where F\ = / )fT /2/ and 

F 2 = ~ff,a/2F. 



Note that the model with real d(o~) satisfy the Weak Energy Condition whereas for 
imaginary d(a), G AB u A u B become a complex quantity. 

2. f( a ) = sechV) 

The geodesies in presence of the decaying warp factor are shown in Fig. [51 As mentioned 
earlier the non-symmetric components drastically change the energy density profile compared 
to the symmetric scenario. So it is expected that the geodesies to be modified. Fig. [5] 
shows that, for d(a) = 0, the geodesies run away from a = and become finite with 
increasing rj (similar results were found in [28]). Whereas in presence of the non-symmetric 
components, geodesies become confined within finite distance from a = 0. This behaviour 
can be understood (in the same line of reasoning as before), by comparing the effects of 
the coefficients appearing in the geodesic equation, from Fig. El Thus in our construction 
of thick branes in a non-symmetric bulk, time-like geodesies are confined for both growing 
and decaying warp factors. As expected in the case of decaying warp factor, with imaginary 
d(a), G AB u A u B is complex. 
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f(^)=sechV). L |,;,r| ± .ScLh-! l r)-2cosh 2 ((r)|~: 



t(cr)=sech V). d(ir)=±j[2cosI) V)-sechV)] 3 



x(A) 
<r(A) 



x(A) 
<r(A) 



fdr) >cdr(<r). d(iT)=0 



>7(A) 



17(A) 



17(A) 



FIG. 5: Plots of geodesies Xj(A) and cr(A) vs t/(A) for decaying warp factor with d(a) real, imaginary 
and zero respectively. Initial conditions are Xj(0) = cr(0) = 0, r](0) = and Xj(0) = d(0) = 1, 7/(0) = 
>/5. 



f(cr)=sechV), d(cr)=± 


5sech 4 (cr)-cosh 2 (o-)J 2 


3 
























-2 
-3 



f(o-)=sech 2 (o-), d(o-)=±iL2cosh 2 ( ( r)-sech 4 (o-)J" 



F 2 



-2 -1 



1 2 



FIG. 6: Comparing the coefficients appearing in the geodesic equation where F\ = /, CT /2/ and 

F 2 = ~ff,a/2F. 



IV. DYNAMIC BULK SPACETIME 

In the previous section, we showed that the non-symmetric components can effectively 
confine the test particles at the location of the brane.In this section, we study the cosmolog- 
ical scenarios in the context of 5 dimensional non-symmetric models. We assume here that, 
a = is the location of the brane and the 5-dimensional spacetime is given by 





(- 


0(77) 








Kv)\ 







a(rj) 








Kv) 


9ab = 








a(rj) 





Kv) 













0(77) 


Kv) 




v- 


Kv) -Kv) 


-Kv) 


-6(77) 


r 2 j 



(4.1) 
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where the 4D part contains the usual cosmological scale factor a(rj) and b(rj) is the corre- 
sponding non-symmetric component along the extra dimension. Here rj is the conformal 
time. The corresponding line element is given by 

ds 2 = a{rj) \-drf + dx 2 } + r 2 do 2 . (4.2) 

Note that, in this model, there is no notion of location of the localised energy density or 
brane. 

Solving Eq. (12. 2p for the above metric (14. ip . the components of as function of a(rj) 
and biff) can be computed (see Appendix). Using these the constraint equations (12.101) lead 
to 

^a v + r 2 (^ + b -A =0. (4.3) 



a 2 * \a b f 

One can easily notice the similarity between Eq. (14. 3 p and Eq. ( 13. 3p . Here again we note 

that the results do not depend on the sign of b{rf). 

The geodesic equations, obeying the timelike constraint, are given by 

d , 9 9n r 2 a 2 bb „ + 2b A a „ — 2ab 3 b r , , v-^ n9 
" + tJ» + ^> ?) + rV-L* - ( " + E Xt) 

i i 

9 r 2 a 2 66„ — 2b A a v — 2r 2 ab 2 a n — 2ab 3 b „ . 9 „ , J jS 

+ r ; 9,, = 0, (4.4) 

a r 2 a 2 bb „ + 26 4 a „ — 2ab 3 b ri , 9 v-^ \ 
X < + Ta" X ' + rV-L* - ( " + " 5» = °' (4 ' 5) 

2^2/ 



r~a 2 b „ + 2b 3 a „ — 2ab 2 b „ 
a 



arj 2 + 2b ^(Y^ 



Xi + a 



0. (4.6) 



2r 4 a 3 - 8a6 4 

Eqs (I4.4p - (l4.6p reduces to standard FLRW scenario when b{rj) = 0. We study the energy 
density and the geodesies for the two cases- A) the de Sitter brane: i.e. a{rf) cc\/rf where 
i] G [—1, 0] and B) power law brane i.e. a(r]) oc r] m where t] e [1, oo). 



A. de Sitter brane 

Substituting for a{rf) = 1/rf (where the proportionality constant is set to unity) in Eq. 
OO]), we get, 

2 

b(rj) = - rV (4.7) 
v/C 3 r 2 - 2 V e 

where C3 is an integration constant. If C3 > 2/r 2 , 6(77) is real for all rj. However, for 
C 3 < 2/r 2 , 6(77) is singular at rj — (C 3 r 2 /2)^. In Fig. we compare the energy densities for 
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cases C3 = 5 and C3 = 0.01 with symmetric spacetime (i.e. b{rj) = 0). The non-symmetric 
property introduces an initial peak in the energy density which is absent in the symmetric 
case. However the asymptotic behaviours are similar in all the situations. The expression 
of energy density contains (r 4 a 2 — 46 4 ) in the denominator which also appear in the Eqs 
(I4.4H4.6I) . This leads to another singularity in the energy density apart from when 6(77) itself 
diverges. 




n v n 



FIG. 7: Plots of energy densities vs 77 for de Sitter brane with 6(77) real (C3 = 5), singular (C3 = 
0.01) at a particular r\ and zero respectively. 



8(17) = l/jf, b(i;) = rfKS-lifY- 



aft) ~ lAf , b(j?)=0 



FIG. 8: Plot of geodesies Xj(A) and cr(A) vs 77(A) for de Sitter brane with b(rj) real and zero 
respectively. Initial conditions are ccj(0) = cr(0) = 0,77(0) = —1 and ±i(0) = <r(0) = 1,77(0) = y/b. 



The corresponding geodesies, shown in Fig. |8]with specific initial conditions, also behave 
similarly for the symmetric and non-symmetric models. Motion along the extra dimension 
slows down as 77 increases whereas the components along X{ increases monotonically. Note 
again that the geodesies are insensitive to initial conditions. 
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B. Radiation dominated brane 

Substituting for 0,(77) = V m (here again, the proportionality constant is set to unity) in 
Eq. flO} leads to 

A/C 4 r 2 ?7 3m - 2 

where C4 is the integration constant. For simplicity we choose m = 2, that corresponds 
to a radiation dominated brane. Figs. [9] and [10] shows the variations of energy density and 
geodesies for different values of C4 and comparison w.r.t. b(r]) = case. At early universe, a 
sharp peak in the energy density makes the cosmological evolution strikingly different from 
the symmetric case. However the geodesies are not significantly different in those two cases. 
Features shown are qualitatively similar to the previous model of a de Sitter brane because 
they both represent expanding 4D universes. The apparent distinction, in the asymptotic 
limit, between Figs. [7] and M arises because de Sitter expansion results in a big rip whereas 
the radiation dominated universe becomes flat at large rj. 



a(f?) = t? 2 , b(??) = 77/(577 6 -2)? 




2 3 4 5 



a(i7) = 7] 2 , b(rj) = ^ 6 /10-2)? 




1.0 1.5 2.0 2.5 3.0 

n 



aft?) ~ if, b(r,)=0 




12 3 4 5 



FIG. 9: Plots of energy densities vs rj for Radiation dominated brane with b{ry) real, singular at a 
particular 77 and zero respectively. 

Note that one can further analyse these solutions for different values of m which corre- 
spond to different cosmological scenarios. 

V. SEARCH FOR WARPED DYNAMIC BULK SPACETIME 

Let us investigate the possibility of constructing thick warped cosmological 3-branes em- 
bedded in 5D non-symmetric spacetime within the framework we have constructed here. We 
choose the following ansatz: 
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afo) = rf, b(;,) = r,K5rf-2fi 

.DP — 1 — 1 — 1 1 — i — 1 1 — 1 — 1 1 — 1 — 1 1 — 1 — i 1 — 1 — 1 1 — i 



afr) ~ g 2 , b(i;)=Q 



FIG. 10: Plot of geodesies Xj(A) and cr(A) vs 77(A) for Radiation dominated brane with b{rj) real and 
zero respectively.Initial conditions are Xi (0) = <t(0) = 0,77(0) = 1 and x;(0) = <j(0) = 1,77(0) = VE. 



( -a(r))f{a) b{q)d(a)\ 

a(r])f(a) Kv)d(a) 

a( V )f(a) 6(^)^W 

a{r])f{a) b{rj)d{a) 

\-b(r))d(a) -b{ri)d{a) -b{r])d{a) -b(rj)d(<r) r 2 



with the line element 

ds 2 = a(rj)f(a)[-dr] 2 + dx 2 } + r 2 da 2 . 
Then the constraint Eqs. (I2.10p leads to the following coupled equations: 



3— a iV + r 
d 2 

7 



d 2 ) \ a b 

f,a d^ a 

7 ~d 



0. 

0. 



Only possible scenario where Eqs ( 15. 3j) and ( 15. Ah can be solved is when 



a(r]) oc b 2 (r]) 



and 



/(<7)ocd». 



Substituting Eq. ([53]) in Eqs (J53J and O leads to: 

a(rj), b(r]), f(u), d(a) = constants, 



(5.1) 



(5.2) 

(5.3) 
(5.4) 



(5.5) 



(5.6) 



i.e. a flat bulk spacetime. Thus thick brane solutions in a dynamic bulk spacetime is difficult 
to find in this construction. However this difficulty arises from the additional assumption we 
made, through Eq. (12.111) . to avoid solving the recursive equation (12.81) . This assumption 
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essentially cuts down the full solution space to one of its subclasses. Thus it is possible that 
the warped dynamic braneworld geometries exist in the full solution space allowed in the 
framework of non-symmetric geometry and can be derived using the same ansatz used here. 
Another interesting possibility would be to assume different functions as non-symmetric 
components. We hope to report on these aspects in a later article. 

VI. DISCUSSION 

Non-symmetric gravity models have been proposed as a natural modification or extension 
of GR. In this work we have analysed, in detail, a yet unexplored aspect of non-symmetric 
geometry where warped thick branes are embedded in 5D spacetime. Here we summarise 
the key results. 

1. We first introduced the field equations corresponding to non-symmetric geometry. 
We have considered models for which the anti-symmetric Ricci tensor (12. 6p vanishes. 
This assumption considerably reduces the mathematical complexity in solving the field 
equations. On the other hand it leads to a single constraint equation for the metric 
functions. To understand the physical implications of these models we have analysed 
the corresponding energy densities and geodesies. 

2. These models provide a natural mechanism for geodesic confinement. They also satisfy 
the Weak Energy Condition when the non-symmetric metric components are real. 

3. We have constructed stationary 5D bulk spacetime where thick branes with growing 
or decaying warp factors are embedded. The non-symmetric components are derived 
from Eqs. (I2.10p w.r.t. chosen warp factors. The timelike geodesies are analysed in 
both the cases. It is found that the geodesic component along the extra dimension is 
confined, by the geometry itself, near the location of the brane in both the scenarios 
with growing and decaying warp factors. This result is different from the corresponding 
symmetric spacetime cases where such confinement can only be achieved for growing 
warp factors. 

4. In the dynamical scenario, non-symmetric components are derived assuming standard 
FLRW brane embedded in a 5D spacetime. We considered de Sitter and radiation 
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dominated expansion in the 4D spacetime. The corresponding energy densities show 
considerable departure from the symmetric models at early times. This suggests that 
the early universe dynamics is different. It may be interesting to look at the scalar 
power spectrum for these models |30|. However the geodesic profiles are not drastically 
modified by the effect of non-symmetricity. 

To understand the role of non-symmetric properties of the spacetime on the curvature we 
need to study geodesic deviations. Further the geodesic and Raychaudhuri equations need 
to be solved simultaneously to investigate the kinematical quantities: expansion, rotation, 



and shear 



311 ] . Such models could be constrained using the observational data |l9j. We 



hope to report on these aspects elsewhere. 
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Appendix 

The affine connections derived from Eq. ( 12.2ft for stationary and dynamic models are 
given below. 

1. Stationary model 



r 



-pO _ -nO 
1 04 ~~ 1 40 



"/- 2r 2 f 

r 2 ff, a + Zd 2 f, a - 2 fold. 



2f(r 2 f + 2d 2 ) 



■nO _ -pO _ p4 _ p4 
1 iA 1 Ai 1 ufi nu 



-r 



i«4 



d{df, a - 2fd a ) 
2/(r 2 / + 2d 2 ) 



v ^ fl, 



■p2 ■p^ 

1 Ai ~ 1 i4 



r 2 ff, a + 3d 2 l a + 2fdd, a 4 / i<y d(df, a - 2fd, a ) 

i 1 oo — - n + 



2f(r 2 f + 2d 2 ) 
4 f(Add !(T + r 2 f :a ) 
a 2r 2 (r 2 f + 2d 2 ) ' 



2r 2 r 2( r 2f + 2( p) 



(6.1) 
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2. Dynamical model 



r° 

1 oo 



1 04 



1 i4 
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V p0 pO pO 

2a' ' ; " '" 



r o 

1 Oi 



a(r 4 a 2 - 46 4 ) 
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1 40 
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1 4i 
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1 2a 2 (r 4 a 2 - 46 4 ) 

r 2 a 2 b(r 2 a jV + 566 ;T7 ) - r 4 a 3 &,„ - 2ab 3 (r 2 a iTI + 366 ;T7 ) + 2& 5 a,„ 

2a 2 (r 4 a 2 - 4o 4 ) : 



^ r 2 6(r 2 a 2 6 jr; - 2o 3 a - 2ab(r 2 a^ + 66 j7? )) 
44 = ! a 2 (r 4 a 2 - 4o 4 ) ' 



p« _ 
1 oo — 



b{r 2 a 2 b^ + 2o 3 a - 2a6 2 6 fl ) 



1 Oi — 1 iO 



1 Oj - 1 jO 



a(r 4 a 2 - 46 4 ) 
r 2 a 2 {r 2 a^ + — 2b 3 (ba^ + afr,,) 
2a(r 4 a 2 - 4o 4 ) : 
b{r 2 a 2 b. n + 2b 3 {ba. q - a&,„)) 
2a(r 4 a 2 - 46 4 ) ' 



-r 



i4 



b(r 2 a 2 (r 2 a tri + 66^) - 2b 3 (ba jV + ab tV )) 
2a 2 (r 4 a 2 - 4o 4 ) ' 



pi 



P« 

1 04 



_ . _ 6 2 (r 2 a 2 6, r? + 26 3 (5a,-a&,)) 

i j4~ r,~2/„4~9 ,(L4\ i L TJi 



-r 



40 



2a 2 (r 4 a 2 - 46 4 ) 
r 4 a 3 o iJ? - r 2 a 2 b(r 2 a tV + 366^) - 26 4 (6a >r , - ao >r? ) 
2a 2 (r 4 a 2 - 46 4 ) 



p4 
1 iO 
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2a 2 (r 4 a 2 - 46 4 ) 
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